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Q\ , Abstract 

The simplest unified extension of the Minimal Supersymmetric Standard Model with 
bilinear R-Parity violation provides a predictive scheme for neutrino masses which 
can account for the observed atmospheric and solar neutrino anomalies. Despite the 
smallness of neutrino masses R-parity violation is observable at present and future 
high-energy colliders, providing an unambiguous cross-check of the model. 

P 



1 Introduction 



The announcement of high statistics atmospheric neutrino data by the SuperKamiokande 
collaboration [1J] has confirmed the deficit of muon neutrinos, especially at small zenith 
angles, opening a new era in neutrino physics. Although there may be alternative solutions 
of the atmospheric neutrino anomaly it is fair to say that the simplest interpretation 
of the data is in terms of to v T flavor oscillations with maximal mixing. This excludes a 
large mixing among v T and v e 0, in agreement also with the CHOOZ reactor data [^]. On 
the other hand the persistent disagreement between solar neutrino data and theoretical 
expectations Q has been a long-standing problem in physics. Recent solar neutrino 
data H are consistent with both vacuum oscillations and MSW conversions. In the latter 
case one can have either the small or the large mixing angle solutions, with the latter 
being clearly preferred |fj. 

Many attempts have appeared in the literature to explain the data. Here we review 
recent results obtained in a model which is a simple extension of the MSSM with 
bilinear R-parity violation (BRpV). This model, despite being a minimal extension of the 
MSSM, can explain the solar and atmospheric neutrino data. Its most attractive feature 
is that it gives definite predictions for accelerator physics for the same range of parameters 
that explain the neutrino data. 

1 To appear in the Proceedings of the Corfu Summer Institute, Corfu 2001, Edited by J. Rizos. 
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2 Bilinear R-Parity Violation (BRpV) 



2.1 The Model 

The superpotential W is given by 

W=e ab [h^U 3 H b u + hfiQtDjH^ h%l\R 3 E a d - ^H a d H b u + ^M] (1) 
while the set of soft supersymmetry breaking terms are 

Vsoft = MfQTQ) + MfU*Uj + M% 2 D*Dj + M$ 2 L?L? + MfR*R 3 
+m 2 H HTH a d + m 2 H H a u *H a u - W \M.XX + e ab (A^U^ 

. i 

+A%Q h fi 3 H« + A%l b t R,H a d - B^H a d H b u + B^H*) + h.c (2) 

The bilinear R-parity violating term cannot be eliminated by superfield redefinition. The 
reason is that the bottom Yukawa coupling, usually neglected, plays a crucial role in 
splitting the soft-breaking parameters B and B t as well as the scalar masses m? H and Ml, 
assumed to be equal at the unification scale. The BRpV model is a 1(3) parameter (s) 
generalization of the MSSM. It can be thought as an effective model showing the more 
important features of the Spontaneous Broken R-parity model (SBRP)|], |TU[ at the weak 
scale. The mass matrices, charged and neutral currents, are similar to the SBRP-model 
if we identify 

e = v R h u (3) 
The model has the MSSM as a limit when e« — > 0. 

2.2 Radiative Breaking 

At Q = Mqut we assume the standard minimal supergravity unifications assumptions, 

A t = A b = A T = A , B = B 2 = A - 1 , 

m 2 Hd = m 2 Hu =M 2 = M 2 R = M 2 =M 2 =M 2 D =m 2 , 

M 3 = M 2 = M 1= Mi/a (4) 

In order to determine the values of the Yukawa couplings and of the soft breaking scalar 
masses at low energies we first run the RGE's from the unification scale Mqut ~ 10 16 
GeV down to the weak scale. For details see || [Tfl. 

3 Tree Level Neutrino Masses and Mixings 
3.1 Neutral fermion mass matrix 

In the basis ip 0T = (— iX', —i\ 3 , H d , H^, u e , u^, u T ) the neutral fermions mass terms in the 
Lagrangian are given by 

C m = - l -(^) T M N ^ + h.c. (5) 



where the neutralino/neutrino mass matrix is 



M 



N 



with 



where a,- 



1 g' v ii \g v ii 0, Cj). This neutralino/neutrino mass matrix is diagonalized by 



Mi 






M 2 



in 



-\g'vd \g'v u 
\gvd -\gv u 



(6) 



-\g'v d \gv d 
\g'v u -\gv u 



o 



7' 







a 3 



(7) 



^MnM- 1 = diag( 



m x o , m x o , m x o , m x o , m Ul , m„ 2 , t 



3.2 Approximate diagonalization of mass matrices 

If the ]fl p parameters are small, then 

£ = m ■ M't < 1 



(9) 



one can find an approximate solution for the mixing matrix Af. Explicit expressions can 
be found in Ref . B . In leading order in £ the mixing matrix M is given by, 



N* W 1 - i£t£ £t 

o k t j \ -e i - ke f 



(10) 



This decomposition block-diagonalizes Mn approximately to the form diag(m e //, M. x o) 
where 



x T M l9 2 + M 2 g 
m e ff = -m-M m - 



Adet(M x o 



t2 / A 2 A e A M A e A T \ 
A e A M A^ A M A r 
V A e A T A M A T A 2 T ) 

The sub matrices N and V u in Eq. (jlO|) diagonalize A4 x o and m e // 

N*M. x oW — diag(m x o) ; Vjm e ffV u = diag(0, 0, m„), 

where 



111 



m„ = Tr(m e ff) 



M l9 2 + M 2 g 
4det(M x o) 



/2 



■|A| 2 . 



(12) 
(13) 



4 One Loop Neutrino Masses and Mixings 
4.1 Definition 

The Self-Energy for the neutralino/neutrino is 
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I < p 



Then 
with 
where 
and 



AM 



I (nJ5K a )+nS(mJ))-| (m x oSK(m?)+m x oEg(m| 



A=0 



K 531 



e*) ; n y = |(n L + u R ) 

2 

: - 7^ + In 47r 

4 — a 



(14) 

(15) 
(16) 
(17) 
(18) 



4.2 Diagrams Contributing 

In a generic way the diagrams contributing are 



P' \ 
i i 



? f %q 

l P J 



p-q 

g :' \ q 
->—> — > — ■-»- 



p-q 



q : \ q 
->— ' — * — <->- 



p j i p j 



These diagrams can be calculated in a straightforward way. For instance the W diagram 
in the £ = 1 gauge gives 



16n 



l — 2 £ 2 (Of^OSJ + O^OS) Brip 2 , ml m 2 w ) 



= ~T^~2 E("4) {OffiOft + 0^02 n ,T) m fe £ (p 2 , ml, n&) 
iD7r fc=i 

where Bq and I?i are the Passarino-Veltman functions, and O cnw , O ncw are coupling 
matrices. Explicit expressions can be found in M . 

4.3 Gauge Invar iance 

When calculating the self-energies the question of gauge invariance arises. We have 
performed a careful calculation in an arbitrary gauge and showed || that the result 
was independent of the gauge parameter £. 
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4.4 The One-Loop Mass Matrix 

The one-loop corrected mass matrix is 



diag 

where 

A, 

diag 

Now we diagonalize the 1-loop mass matrix 



M 1L = M ^ + AM 1L (19) 



MdL = MM N M T (20) 



Mdiag — N'M 1L N' T (21) 

Then the mass eigenstates are related to the weak basis states by 

mass » weak ( 99\ 

AO ~~ Jv ia AO 

with 

M lL = M' M (23) 

The usual convention in neutrino physics 

v a = U ak v k (24) 

is recovered in our notation as 

U ak = A/" 4+M+Q , (25) 

4.5 Approximate Formulas for Masses and Mixings 
4.5.1 The masses 

Looking at the numerical results |J we found that the most important contribution came 
from the bottom-sbottom loop. To gain an analytical understanding of the results we 
expanded the exact results in the small ftp parameters. The result is [11 

/ AiAi A X A 2 AiA 3 \ I exEi e x e 2 e x e z 



where 



A 2 A X A 2 A 2 A 2 A 3 + ci e 2 e 1 e 2 e 2 e 2 e 3 | (26) 
VA 3 Ai A 3 A 2 A3A3/ \e 3 e! e 3 e 2 e 3 e 3 



,2 \ 1 ,,'-2 > 1,2 m 2 - 

4det(M x o) ' " l 16tt 2 > 2 "°^ 



M l9 ' + M 2 g» 3 . h . - b2 

c o = , , 7777 x ; ci = — — m fe sin 26*6 — log —f (27) 



Diagonalization of the mass matrix gives [11 



(28) 



3 ™ 6 sm2,4log!l^ (29) 



U2 ~ 16tt 2 V rn\ |A| 2 
m - * 4det(M x o) |A| (30) 
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The formula for m„ 3 is the tree-level formula that we used to fix the scale of the at- 
mospheric neutrinos by choosing |A|. The quality of the approximation formula for m U2 
can be seen in Fig. |I[ Only for small tan (3 and very small b, b mixing the approxima- 
tion can lead, in some cases, to incorrect results. Details of the derivation can be found 
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Figure 1: a) Ratio of m^ p , given by Eq. (|29|), with the exact numerical result, b) The same 
with the cut sin26 b > 0.02. 



in Ref. |TT[ where the second most important contribution, coming from the loop with 
charged Higgs/charged leptons, is also discussed. 



4.5.2 The mixings 

The atmospheric angle is easily obtained in terms of the ratio A 2 /A 3 . For the solar angle 



in the same approximation we also get a simple formula [[1 1 

2e 2 



tan 2 9. 



ol 



(62 + 6 3 ) 2 

that is also in very good agreement with the exact result. 



(31) 



5 Results for the Solar and Atmospheric Neutrinos 

5.1 The masses 

The BRpV model produces a hierarchical mass spectrum for almost all choices of param- 
eters. The largest mass can be estimated by the tree level value using Eq. fl3"0p. Correct 
Am 2 tm can be easily obtained by an appropriate choice of |A|. The mass scale for the solar 
neutrinos is generated at 1-loop level and therefore depends in a complicated way in the 
model parameters. However, in most cases the result of Eq. ([29|) is a good approximation 
and there is no problem in having both Am 2 tm and Am 2 oiar set to the correct scales. 

5.2 The mixings 

Now we turn to the discussion of the mixing angles. We have found that if e 2 /|A| 
100 then the 1-loop corrections are not larger than the tree level results and the flavor 
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composition of the 3rd mass eigenstate is approximately given by 



U, 



a3 



Aa/|A| 



(32) 



As the atmospheric and reactor neutrino data tell us that — > v T oscillations are preferred 
over — > v> e , we conclude that 

A e < A M ~ A r (33) 

are required for BRpV to fit the data. This is sown in Fig. a). We cannot get so 
easily maximal mixing for solar neutrinos, because in this case U e3 would be too large 
contradicting the CHOOZ result as shown in Fig. ^| b). 
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Figure 2: a) Atmospheric angle as a function of lA^I/yAj + A£. b) as a function of 
\A e \/jAl + A|. 



We have then two scenarios. In the first one, that we call the mSUGRA case, we have 
universal boundary conditions of the soft SUSY breaking terms. In this case we can 
show || that 

7 = T < 34 > 

Then from Fig. ^| b) and the CHOOZ constraint on U^ 3 , we conclude that both ratios in 
Eq. ( |3"1D have to be small. Then from Fig. |3| we conclude that the only possibility is the 
small angle mixing solution for the solar neutrino problem. In the second scenario, which 
we call the MSSM case, we consider non-universal boundary conditions of the soft SUSY 
breaking terms. We have shown that even a very small deviation from universality of the 
soft parameters at the GUT scale relaxes this constraint. In this case 



£i -L — 

£„. A„ 



(35) 



Then we can have at the same time small U^ 3 determined by A e / A^ as in Fig. ^| b) and 
large tan 2 (6' so ; ar .) determined by e e /e M as in Fig. ||b). 
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Figure 3: Solar angle as function of: a) |A e |/\ k 2 + A 2 . ; b) e e /e 



6 Probing Neutrino Mixing via Neutralino Decays 

If R-parity is broken, the neutralino is unstable and it will decay through the following 

Itqq', 



channels: Xi 



ViVjV kl Uiqq, 



v, If I 



j L k ) 



Hi- 



lt was shown in Ref. [12 1, that the 



neutralino decays well inside the detectors and that the visible decay channels are quite 
large. This was fully discussed in Ref. [12| and is illustrated in Fig. |]. We have seen 
before that the ratios |Aj/Aj| and | / e^- 1 were very important in the choice of solutions for 
the neutrino mixing angles. What is exciting now, is that these ratios can be measured in 
accelerator experiments. In Fig. |5| a) we show the ratio of branching ratios for semileptonic 
LSP decays into muons and taus: BR(x n'o)/BR{x — * Tq'q) as function of tan 2 9 atm . 
We can see that there is a strong correlation. In Fig. [5] b) is shown the ratio of branching 
ratios for semileptonic LSP decays into muons and taus: BR(x eq'q)/BR(x — > ^q'q) 
as function of U^ 3 . Again we obtain a strong correlation. The spread on those figures 
can in fact be explained by the fact that we do not know the SUSY parameters. This 
is illustrated in Fig. || where we considered that SUSY was already discovered with the 
following values for the parameters, 



M 2 = 120 GeV, n = 500 GeV, tan = 5, m = 500 GeV, A = -500 GeV 



(36) 



7 Conclusions 

The Bilinear R-Parity Violation Model is a simple extension of the MSSM that leads to a 
very rich phenomenology. We have shown that the radiative breaking of the Gauge Sym- 
metry can be achieved in this BRpV Model. We have calculated the one-loop corrected 
masses and mixings for the neutrinos in a completely consistent way, including the RG 
equations and correctly minimizing the potential. We have shown that it is possible to 
get easily maximal mixing for the atmospheric neutrinos and both small and large angle 
MSW. We emphasize that the lightest neutralino decays inside the detectors, thus lead- 
ing to a very different phenomenology than the MSSM. If the model is to explain solar 
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and atmospheric neutrino problems many signals will arise at future colliders. These will 
probe the neutrino mixing parameters. Thus the model is easily falsifiable! 
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